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Charged Black holes in Gauss-Bonnet extended gravity are studied. The electromagnetic 
field is coupled non-minimally, as in U(2, 2) Chcrn-Simons theory. We find that the geomet- 
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rical properties of the solution exhibit "phase transitions" as one varies the mass and charge. 



The full phase diagram for all values of the ADM mass and charge is displayed. 
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<N ; I. INTRODUCTION, MAIN RESULTS AND CONCLUSIONS 

In a five dimensional Universe, the Gauss-Bonnet density y/—g(R 2 — 4:R^ U R^ v + R^ vXp R^ u \ p ) 
cannot be omitted in the gravitational Lagrangian. This term is covariant, its associated Einstein 
tensor is conserved, and, despite being quadratic in the curvature tensor, yields second order field 
equations for the metric [J, y] • The most general action for Gravity in five dimensions is then 

= / v/=? ao + ai R + a 2 (R 2 - iR""R^ + R^'R^xp) (1) 

Q_|' JM 5 1 J 

The presence of this term of course changes the dynamical equations, and many aspects of general 
relativity have to revisited. This issue is particularly relevant in the context of brane worlds models, 
$—1 ' and many papers have recently been devoted to the subject 



The simplest problem that can be analyzed in a closed form is the spherically symmetric five- 
dimensional black hole spacetime, 



ds 2 = - N 2 d f + ^l + r 2 dn 3 . (2) 



P 

Although the equations can be solved for arbitrary values of the three couplings 



4. 



we are 



interested in the rol of the Gauss-Bonnet term and then we set, for simplicity, a\ = 1, ao = and 



«2 > 0. The equations of motion yield 



IB Q 



N 2 



f = 1 + — (r 2 - yV + Aa 2 M] (3) 

(To V / 
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where M is an integration constant, that will be seen to be the ADM mass. 
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We first note that for M = the metric reduces to flat space, which is the stable background 
for this theory ^j. It also follows that for r 4 3> AMot2, 



(4) 



showing that, asymptotically, © approaches the five dimensional Schwarzschild metric. Inciden- 
tally, note that if M > «2 then r 4 3> 4Ma2 would hold all the way to the Schwarzschild horizon 
at r + = y2M. If «2 = 0, © becomes an exact solution. 

The function (J3J) has some interesting properties. The value of iV 2 at the singularity r = is 
finite, 



V «2 

although the curvature is still singular 16]. The location of the horizon, N 2 (r+) = 0, is 



r+ = V 2M "T- (5) 

We see from this expression that the horizon exists only for M > a2/4. We thus find a mass gap 
separating flat space from the spectrum of black holes: 

M = Flat space 
< M < a 2 /4 Naked Singularities (6) 
M > a 2 /4 Black holes 

The mass gap appears in all odd dimensional theories containing the highest Lovelock Q] term. In 
three dimensions, this term is just the Hilbert term and the mass gap is present ^J. In this case, 
however, the "naked singularities" have a sensible interpretation in terms of particles ^]. 

In this paper we add electric charge q to this black hole and study the corresponding spectrum. 
We shall see that the solution has some peculiarities not present in usual charged black holes. 

As a first surprise, the addition of charge does not imply the existence of two horizons. There 
are open regions in the plane {M, q} having non-extremal black holes with one horizon. In order 
to find solutions with two horizons, q has to be bigger than a certain critical value, q > q c . Fig. ^ 
gives a summary of the properties of various solutions obtained by varying the values of M and q. 

In one region of the space of parameters, we find a set of non-extremal charged black holes 
having only one horizon. We call this region the "Heavy Branch" because it is defined by the 
condition M > M cr u. For masses within the range 

M crit >M> M exU q>q c (7) 



3 



Muss Gup 
Hal Spaa 



The "Phase Diagram" 



Heavy Branch 

ONE HORIZON 



q=q c 

Critical Charge 



Critical Curve 

M = M cri , 
\ 




Exlreme black holes 
M = Mext 



NAKED SINGULARITIES 



FIG. 1: The phase diagram. For each value of q > q c there exists two black holes phases. Flat space is 
located at q = M = 0, and it is disconnected from the black hole spectrum by a set of naked singularities. 



we find the "Light Branch" with black holes with two horizons. At M = M ex t, we find extremal 
black holes with only one horizon (and zero Hawking temperature). Below the extreme value, the 
solution represents a naked singularity. There is also a "critical charge" below which the light 
branch ceases to exist. Flat space is located at M = q = 0, and thus the mass gap persists in the 
charged solution. 

The terms "phase structure" and "phase transitions" are used here only in analogy with the 
statistical mechanics concept, without implying a direct connection. Of course, given the thermo- 
dynamical properties of black holes, this may turn out to be more than an analogy but we shall 
not study this issue here. 

To avoid future confusions, we stress that the action considered in this paper is not the usual 
minimally coupled Einstein-Maxwell system (plus a Gauss-Bonnet term). Those solutions were 
studied in Ml and do not exhibit this phase structure. Instead, we consider a five-dimensional 

Q 

Chern-Simons theory for the group U(2, 2) |9(, which has a sensible interpretation as a gravitational 
plus electromagnetism theory, with a Gauss-Bonnet term. This interpretation, however, requires 
a symmetry breaking term because otherwise the equations of motion differ from the usual ones 
even asymptotically. This point was discussed in detail in |lf| . 

The application of Chern-Simons theories to gravity has been discussed several times in the 
literature and we shall not repeat it here. The first constructions were reported in three dimensions 
in and the same idea was then applied in to five dimensions. See Q for other aspects. 

For the purposes of this paper, we refer the reader to Ref. ^(| were many details omitted here 
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can be found. In particular, the asymptotic form of the charged black holes was already reported 
in that Reference. The goal of the present paper is to display the exact solution for an arbitrary 
mass M and electric charge q, and study the associated phase space. 



II. THE EQUATIONS AND THEIR SOLUTION 
A. The equations and spherically symmetric ansatz 

We start by writing down the equations of motion associated to a Chern-Simons theory for 
the group £7(2,2). This group contains SO(4, 2) which would give pure gravity with a Gauss- 
Bonnet term. The extension to {7(2,2) incorporates an Abelian one form, that we interpreted as 
an electomagnetic field, coupled non-minimally to gravity. See for a detailed analysis, and 
other motivations to study this system. 

A key issue in the analysis (and discussed in [if] ]) is the fact that the Chern-Simons equations of 
motion do not provide a sensible theory for the £{7(2, 2) field. For example, there is no linearized 
theory and the metric is not asymptotically Schwarzschild. Given the strong topological roots of 
the Chern-Simons construction, it is still attractive as a field theory, and one would like to know 
if the equations can be "repaired" by some mechanism, hopefully within the same theory. Some 
progress in this direction was reported in 

Here we follow |lfl| in which a symmetry breaking term is added to the action. The Chern- 
Simons equations then becomes closer to the real world and one can start asking questions such 
as what is the structure of black holes, and what is the nature of the couplings between the 
gravitational and internal U(N) gauge fields degrees of freedom. In this paper we concentrate on 
the coupling between the gravitational and {7(1) field. 

The symmetry breaking term added in [10| consists in a cosmological term (vacuum energy) 
and it is parameterised by a real number r. Let e a be the five dimensional vielbein one-form, w ab 
the spin connection, and A the Abelian one-form. The equations of motion following from the 
{7(2,2) Chern-Simons theory (including the symmetry breaking term proportional to r) are, 

e abcde [{R ab + e a e b )4R cd + e a e h )-T 2 e a ^e\e c ^e d ] = -AT £ aF (8) 

e abcde (R ab + e a e b ),T c = (R de + e d e e )*F (9) 
^R ab ,R ab - d(e a ,T a ) = FaF (10) 

where R ab = dw ab + w a c w cb ,T a = de a + w a h e b and F = dA. For r = 0, these equations are 
equivalent to T A T = with T G {7(2, 2) which are the exact Chern-Simons equations. 
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Since we are interested in black hole solutions, we write the ansatz for the metric and gauge 
field with spherical symmetry 

ds 2 = -N(r) 2 dt 2 + f 2 (r)dr 2 + r 2 cM 3 , (11) 
A = -(/){r)dt, (12) 

where N, f and (p are functions to be determined. 

As shown in the dynamics of the Abelian form A is linked to the torsion tensor. If we 
assume from the very beginning that T a = 0, then the equations of motion for A do not give 
Maxwell's equations in any limit. On the contrary, letting T a to be different from zero yields a 
system of equations that can be analyzed perturbatively and one obtains, to first order, Maxwell's 
theory for the potential (fi(r). The relationship between the torsion T a and the Maxwell field F is 
encoded in Eq. Q. To first order we ignore the right hand side of JHJ) and find the gravitational 
background (AdS space) R ab = (r — \)e a e b . Replacing in © it implies Q 

e a f\T a = *F, (13) 

where * represents Hodge's dual. Replacing (j!3j) into (|10|) one obtains the usual five-dimensional 
Maxwell-Chern-Simons system dF = F A F. This is, in short, the mechanism that transforms the 
"topological" one-form A into a radiating Maxwell field. 

Incorporating the back reaction from the right hand side of (jHJ) generates corrections to (|13[1. 
and the Maxwell equation. It is precisely the role of these corrections what we aim to investigate 
in this paper. 

Let us then assume that the torsion is not zero, and let w ab = w ab (e) + K ab where w ab (e) is the 
solution to the equation de a + w a b (e)/\e b = 0, and only depends on the metric, k parameterizes a 
non zero torsion by T a = K a b Ae b . In order to prescribe a spherically symmetric ansatz for K ab , it is 
convenient to express all indices in the coordinate basis n a p ^ = e aa K ab tl e^p ■ The correct ansatz 
with spherical symmetry for this tensor follows by studyin g th e equations of motion order by order 
starting from the AdS vacuum. The details were given in jlOfl. and the result is, 



KflU A = 27K Va Q/? z <*0 ~ z ^ Ux + 2 9X[ixV u ] (14) 
where e ^ AQ/3 is the Levi-Civita tensor with e * r010203 = 1, \fh = r 3 sin 2 6 1 sm9 2 , and z = dUdr, U = 
P(r) dt, V = a{r) dr. The torsion then contributes with three more functions of r to be determined 
by the equations of motion, namely, a(r),f3(r) and ip( r )- ^ is closely related to the electric field 
(see Eq. I)13|l). while a and (3 are auxiliary fields which will be eliminated algebraically from their 
own equations of motion. 
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We shall see that this ansatz does provide an exact solution to the full U(2, 2) system. 



B. The spherically symmetric reduced equations 

In this section we present some of the details in finding the equations and their solution. The 
reader interested only the final result can jump to the next section where the final form of the 
metric is displayed and its properties analyzed. 

It is a direct but long calculation to replace in the equations of motion 1)8110(1 the ansatz with 
spherical symmetry shown in the last paragraph. The resulting equations turns out to be extremely 
complicated. Some simplification can be achieved by making field redefinitions which simplify the 
expressions for the curvature R ab and torsion T a . These field redefinitions involve all variables. 
We transform {a,/3, N,f, ip,4>,} — > {oti, /3i, N%, f2, i^i^i^ }, according to 

ot\ — 1 
a = 

r 

(3 = N +N' + Na 



Jl 

N = f x N u N 1 = ef N * 

dr (i\qV>i*i) (15) 



tb = — 
r 

Note that the horizon structure will be controlled by the zero'es of the function f\. 

Inserting the above ansatz into the equations of motion, we find two sets of equations which 
can be analyzed and solved separately. 



1. The a,x , ipi , /2 system 

The first set of equations involves only the functions {a\,ipi, /2}. The equations are the following 
(prime indicates radial derivative): 

(#? - 4r 2 ^i - 4^1 - 1M\Wx ~ 2/ 2 r + 2rf 2 a 1 + 2r 2 r 3 ai = (16) 

aih ~ h - 2nM4 = (17) 
- (Vi/ 2 )' = (18) 
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We note that (|T7j) and (|18j) can be easily solved. We find, respectively, 

«, = i+^ (19) 

00) 

where go is an integration constant that will be related to the electric charge. It will be convenient 
in what follows to do yet one more redefinition, 

Mr) = ~^~y qo = — (21) 

where ^(r) is a new function of r. Replacing ()20j) into (|19|) we find the following expression for a±, 

(Pr)' 

ai = ^- (22) 



where P is a short hand for 

\J/3 _ 2 



^3 • 

Finally, we replace (|19|) and (|2U|) in 1)16(1 and obtain a closed equation for the function 



(23) 



(V* 2 + 3g 2 * 6 - 6g 2 r% 4 - 6q 2 * 5 - 4 g 6 + ^- + 2r 2 r 3 tt V — <f 3 ) = (24) 
V / dr 

This last equation is linear in VP'. The integral can be done explicitly and we find the implicit 
solution for ^(r), 



T 



i f^3_ q 2\ 4 1 q 2 4 q 6 g q 6 g q 8 2 4 q 



r 4 \J/ 2 _ 3 _H u -I _H 1. _H u _ -2 1 : -2. L 

2 ^ * 3 J 2 $ 2 7 2 $ 8 10 $ 10 11 * n 

,24 ?1_12^ , 1^1_3 
5 ^ 5 13 4 2 ^ 4 7 ^ 14 ^ ' 

where M is an integration constant that will be seen to be the ADM mass of the solution. This is 
an algebraic equation that should be inverted to find \P(r). 



2. The /3x,N 2 ,$ system 

We now proceed to find expressions for /3i,iV~2 and $i in terms of The three remaining 
equations are: 

= f 2 r + Pif2 + r 3 (26) 
= anh+fife+ftfa + Wife ~ ZM'iPi + 

+2r 2 + 4r^i 2 ^ -2ri\) 1 ^' 1 - 2a 1 j3 1 r + r f 2 + ai r 2 + 2r(3 1 (27) 



= 2 ai Pt - 2fo$t - 2N 2 i> 1 2 + 2N 2 +f 2 + 



2 

2 



+2 N 2 f 2 ~ 2 tpi + 2 r# + 2 r flJ + 2 r + 2 rN 2 Pi + 2 r 2 JVg (28) 
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This set of equations can also be solved in a closed form. We shall not go into the details on 
how to find the solution, we only quote the result. Using the value of a\ found in the previous 
paragraph, we find that Eqns. ()26I28|) are solved by: 

h 

N 2 = ^ (30) 

J2 

1 (3/| - 2r 4 U' , x 

Note that making the redefinition (|21|) . and using the solution (|2()|) for f 2 we have 



N 1 = ef N * 



\J/3 



-,2 



4 



^3 



This completely solve the problem. All functions are known in terms of r and and ^ is known 
in terms of r by (|25|) . 

III. CHARGED BLACK HOLES AND A PHASE TRANSITION 

A. The metric 

Let us summarize the results of the analysis of the equations of motion. The metric ansatz was 

ds 2 = -N(rfdt 2 + f(r)dr 2 + r 2 dn 3 . (32) 
The functions N and / are fixed by the equations of motion as (prime denotes radial derivative) 



N 2 = P 8 fl 



f 2 



(Pr)' 



4 (33) 
Ji 



where 



/i 2 = l + - 2 -^-^ (34) 

and P is given in (|23|) . Finally, ^ is a function of r defined by the algebraic equation (|25|) . 

In some applications it may be convenient to define a new radial coordinate, as suggested by 
Eqns. (t23|) and (IHHl). 

p = Pr (35) 
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In terms of this new coordinate the metric takes a simple form, 



ds 2 = -P 8 f 2 dt 2 + ^ + P'dnfj (36) 



and the relation (|23|) becomes 

2 4 t2 Q 2 8 q 6 q 6 9 q 8 48 q 8 
r 2 o 4 + AM = * 2 + 6 - h 9 - 

48 q 4 24 q 10 1 q 12 q A 6 g 10 

IT W + 13 ~ 2 ^6 + 3 vJ/4 + 7 xp4 v 37 ) 

From now we shall only consider the case r = 1. This is only for simplicity in some calculations, 
but it does not affect the main conclusions. 



B. Known limits of the solution. Reissner-Nordstrom Solution 

Since the exact solution displayed in the previous paragraph is rather complicated, as a first 
check we analyze how this solution reduces to the known ones, in various limits. We first study the 
uncharged solution, found in jj] and discussed in the introduction. Then we show how in the limit 
of small charges and large radial coordinate r, we recover the usual Reissner-Nordstrom spacetime. 

Consider first the uncharged solution with q = 0. In this case, P = 1 (hence r = p), we can 
solve ^ explicitly, \& = \/r 4 + AM, and obtain for f\ the closed expression, 



f 2 = 1 + r 2 - Vr 4 + AM, (38) 

J], described in the introduction, with 0,2 = 1. 
Consider now the charge as a small parameter, and seek for a perturbative solution in q 2 to the 
equation (|23|) . Let 



representing the uncharged solution 



*(r) = y/r 4 + AM + q 2 hi(r) (39) 



and solve (|25l) keeping only the linear terms in q 2 . We find for the first order perturbation 

2r 4 + 3vV 4 + AM 
(r 4 + AM) 



Hr) = Zi ■ „ J2 ( 4 °) 



Replacing in (|33(1 . and taking the limit r 4 3> AM, we find 

coinciding exactly with the Reissner-Nordstrom spacetime. We can also see that the parameter 
that we called q is in fact the electric charge, up to a normalization. 



10 



C. The origin, the Kasner singularity, and physical range of radial coordinate 

The metric (|3fi|) has a curvature singularity at the origin p = 0. There is another singularity at 
the point where P vanishes (P was defined in (|23jl ) , 

- q 2 = 0. (42) 

At this singularity, proper times are shrink to zero while spacelike separations are stretch to infinity. 
The volume element, however, remains finite, detg ~ P 8 x (P~ 2 ) 4 = 1. For this reason, we call 
this point the "Kasner singularity". 

Both singularities can be shown to be physical in the sense that the components of the curvature 
tensor in an inertial frame diverge. It is then important to ask whether these singularities are 
protected by horizons. 

The analysis of existence of horizons is greatly simplified by noticing that the function \& can 
be used as a radial coordinate. Although the algebraic relation (|37j) between p 4 and Vl/ is quite 
untractable and attempts to invert it explicitly are hopeless, we can in fact show that in the domain 
of interest, it is an invertible function. 

We first note that the derivative of (|37j) can be factorized in the form, 

dp A (I' 3 + 2q 2 - y/6q%)(V 3 + 2q 2 + v / 6^)(^ 3 - q 2 ) 4 

~ w • (43) 

Since ^ 3 = q 2 is a curvature singularity, we do not need to worry about non-invertibility at 
that point. We need to focus on either \E' 3 > q 2 or \l/ 3 < q 2 . Since, asymptotically, \P ~ r 2 is a 
large positive number, the physical domain of the function is 

q 2/3 < ^ < +oo, (44) 

and we explore invertibility of (|37[) on this domain. 

All factors in (|43|) are positive definite in the physical domain except for the first one. We 
would then like to know if the solutions ^ c to the equation ^ c 3 + 2q 2 — \^6q^ c = lie in the 
physical range of the function or not. Let us first note this equation has positive solutions only 
if q < q c := (2/3) 3/2 . Thus, our first conclusion is that the relation 1)37(1 is invertible for q > q c . 

Recall now that the metric has two singularities, ^ = q 2 , and p = 0. Let us call the 
particular value of such that ^(^o) = 0. It turns out (this is most easily done by a graphical 
analysis) that if q < q c the three numbers an d q 2 ^ are ordered according to 
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This means that the non-invertible point ^ c is beyond the origin p = and thus, it does not affect 
the physical domain. 

D. Horizon structure 

Given the form (|36j) of the metric it is clear that horizons will arise whenever the function 
vanishes. As mentioned before, the relation between p and ^ is invertible and we can study f\ as 
a function of ^l. We write here the explicit form of /f in terms of the ADM mass M and the 
charge q, 

\J/6 



i 



i 

2 ^ 2 



(tf3 _ q 2f 



(46) 



/ 1 q 12 6q2 24q™_ 

~ 2 4> 16 + ^ 2 + 13 * 13 ~ 5 ~ 11 ^ n ~ 7^ 7 



9 q 8 48 q 8 



8 q 6 6 g 10 9g 6 3g 4 



48 

y ^5 



4M 



In Fig. |2]we have plotted fi(^) 2 in the domain g 2 /3 < \j/ < DO) f or q = Q } and five different values 
of the mass M. (The picture is actually generic for all values q > q c .) Let us analyze each curve 
separately. 




FIG. 2: The function ff ( , 3>) in the domain q 2 / 3 < ^ < oo for q — 6 and five different values of M 



1. Light black holes, M < M cr n, and extreme black holes 

The lightest case, corresponding to M ex t = 7.057..., represents the extreme black hole. It touches 
the horizontal line once, and its derivative is zero there too. The masses M ex t(q) are defined by 
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the equations 

fl = 0, § = 0. (47) 

These equations can be solve numerically and we have found a linear relation 

M ext (q) ~ -0.30914... + (1.2247...) q. (48) 

(The linear approximation is better for charges q > 2.) This result is remarkable because the 
equations (|47|) defining M ex t(q) form a extremely non-linear system. Note also that the slope of 
the curve approaches 1.2247... = \/3/2, which is precisely the value obtained by the asymptotic 
solution (valid for q/r 2 <C 1), 

2 _ 2M 2g 2 / /2 g V /3 



A word of caution is in order here. This analysis does not imply that the function M ex t{q) is 
exactly linear, for all values of q. We only claim that the linear relation is a good approximation 
for that curve. 

Let us now we move to the curve M2 = 7.8. This looks very much like a standard charged black 
hole, f 2 intersects the horizontal line twice, and thus there are two horizons. 

The black holes discussed so far have one or two horizons, and f 2 diverges as one approaches 
— ► q 2 ^ 3 ; in these cases, the Kasner singularity at ip 3 = q 2 is met before the origin p = 0. 

2. The critical mass M = M cr i t 

If we carry on making the black hole heavier, we reach the curve (for q = 6) M cr n = 8.113... 
where something new happens (we give a close expression for M cr n below). This curve intersects 
the horizontal line only once, and thus it has only one horizon. Also, at the origin, f 2 has a finite 
value. 

To have a better understanding of this case, consider the function f 2 displayed in l|46j). For 
generic values of M there is a explicit singularity at ^ 3 = q 2 . However, if M is fine-tuned such 
that the numerator (square root) vanishes at that point, the pole is cancelled. In fact, one observes 
that the zero in the square root is stronger than the zero in the denominator, and that whole term 
vanishes at ^ 3 = q 2 . The value of f\ at that point is then, 

/?U s ,=i-^ 2/3 > m 
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which is in fact finite. Since M enters linearly in the square root, the value of M, called M cr it, 
such that the square root vanishes at Vl/ 3 = q 2 can be calculated directly, 

M crit = g-^y q 2 ' 3 (33 q 2 / 3 + 26) , (50) 

(13!! = 13 * 11 * 9 • • • 3 * 1), 

Finally, recall that the square root is nothing but p 2 (see Eq. (|37j) ). This means that, by 
definition, at the critical mass M cr n, the origin p = and the Kasner singularity Vl/ 3 = q 2 coincide. 



3. Heavy black holes: M > M cr n 

Let us now increase the value of M above M cr it. We find the curves M4 = 8.5 and M5 = 10. 
These curves intersect the horizontal line only once. The associated black holes then have only one 
horizon, despite being charged. 

In this class of solutions (with M > M cr it) the origin p = is met before the Kasner singularity. 
This is the reason that the curve stops before reaching Vl/ 3 = q 2 . At p = 0, f 2 has a finite value 
(just like the uncharged black hole discussed in the introduction). 



4- The critical charge 

We have seen that the spectrum of black holes is separated into two branches, the heavy branch 
with M > M cr n an d the light branch with M < M cr n. The interphase is defined by the critical 
curve M = M cr a displayed in (|50|) which depends on the charge q. We shall now see that there 
exists a particular value of q, namely, 

2 x 3/2 

3) ■ < 51 > 

for which the light branch produces only naked singularities. 

In fact, going back to Eq. (fl9|) we note that for q = q c , the value of f 2 at ^ 3 = q 2 is zero. We 
have plot in Fig. [31 the function f 2 for q = q c and three different masses, M < M cr n, M = M cr u 
and M > M cr a. We observe that the curve M < M cr n cannot intersect the horizontal line. The 
light branch thus gives rise only to naked singularities. 

For charges q < q c the critical curve is pulled upwards. The light branch will carry on producing 
naked singularities, while the heavy branch will have both, black holes and naked singularities 
depending on M. 
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FIG. 3: The function f\ as a function of * for q = (2/3) 3 / 2 . 



These results are summarized in the "phase diagram" displayed in the introduction, Fig. ^ 
showing the various black holes types for all values of M and q. The most important aspect of that 
diagram, and of this paper, is the existence of two types of black holes for charges q > q c , which 
are continuously connected by varying the parameters M and q. 



So far we have only analyzed the properties of the metric. To compute the value of the elec- 
trostatic potential we first go back to the redefinitions (|15j) . and recall that the functions N\,ipi 
and are known in terms of ^, which is algebraically related to the radial coordinate. The full 
expression for the potential is not very illuminating so we don't display it here. We only quote the 
result in the asymptotic limit r — > oo 



showing as stated above that the electromagnetic potential is asymptotically controlled by Maxwell 
equations. 



E. The Coulomb potential 




(52) 
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